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Abstract 

We present a detailed study of the heating effects in dielectric measurements carried out on 
a liquid. Such effects come from the dissipation of the electric power in the liquid and give a 
contribution to the nonlinear third harmonics susceptibility xs which depends on the frequency 
and temperature. This study is used to evaluate a possible 'spurious' contribution to the recently 
measured nonlinear susceptibility of an archetypical glassforming liquid (Glycerol) . Those mea- 
surements have been shown to give a direct evaluation of the number of dynamically correlated 
molecules temperature dependence close to the glass transition temperature Tg k. 190K (Crauste- 
Thibierge et al, Phys. Rev. Lett 104,165703(2010)). We show that the heating contribution is 
totally negligible (i) below 204K at any frequency; (ii) for any temperature at the frequency where 
the third harmonics response X3 is maximum. Besides, this heating contribution does not scale as 
a function of f /fa, with fcx{T) the relaxation frequency of the liquid. In the high frequency range, 
when f / fa > 1, we find that the heating contribution is damped because the dipoles cannot follow 
instantaneously the temperature modulation due to the heating phenomenon. An estimate of the 
magnitude of this damping is given. 
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I. INTRODUCTION 



Upon fast enough cooling, most liquids avoid cristallization and enter a supercooled state. 
The latter is characterized by the "viscous slowing down" phenomenon jlj , an extremely fast 
increase of the viscosity when the temperature T is decreased. Below the glass transition 
temperature Tg the viscosity is so high that the system is, in practice, a solid. The glass 
transition is a longstanding issue of condensed matter physics, since no structural signature 
has ever been detected around Tg, e.g. the static neutron spectra do not change at Tg, 
contrary to what happens for the cristallization [l| . Significant progresses in understanding 
the physics of the glass transition were made in the last fifteen years, when the heterogeneous 
nature of the dynamics of supercooled liquids was estabhshed, through various experimental 
{2-5] as well as numerical works The important concept of 'dynamical heterogeneities' 
has emerged , according to which the relaxation comes from the collective motion of 
groups of Ncorr molcculcs. All these groups evolve in time, some of them being faster or 
slower than the average dynamics. The viscous slowing down would come from the fact 
that this .umber iV„„ ..creases as T decreases towards T,. The temperature dependence 
of Ncorr thus bccamc a crucial issue in the field, triggering new theoretical ideas |7H10|. 
One of them came from an analogy 10| with the well known spin glass physics in which 
the nonlinear suceptibility diverges at the critical temperature T^, reflecting the long range 
amorphous order which sets in at T^. It was proposed in Ref. 10|, that a similar effect occurs 
in supercooled liquids, with the key difference that the peak of the nonlinear response Xsi^) 
should appear at finite frequencies u of the order of l/r(T) where t(T) is the relaxation time 
of the supercooled liquid. For the first time, a (nonlinear) susceptibility was directly related 
to the quantity Ncorr of interest, allowing to scrutinize its temperature dependence. We have 
performed recently the corresponding experiment for the nonlinear dielectric susceptibility 



of glycerol 
decreases. 
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12| and shown that the Ncorr{T) dependence was indeed an increase as T 



The nonlinear susceptibility X3 measured in Refs. [Uj, ll2|] corresponds to the detection 
of the third harmonics of the polarisation P3, at three times the frequency of the applied 
electric field Ecos{u)t). For ~ 1 MV/m, one finds that P3/P1 is typically around 10~^, 
with Pi the polarisation at the frequency u (mainly dominated by the linear response of the 
system). Due to the small value of P3/P1, a thorough analysis is required to see whether 
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some 'spurious' effects can affect the P3 measurements. Expanding on the arguments given 



in Ref. 



13| . we study in this paper the contribution of the ac heating of the hquid to our 



P3 measurements. This heating comes from the fact that the strong apphed electric field 
leads to some dissipated electrical power p with a d.c. component as well as a component 
P2{t) oscillating at 2u). Before being absorbed by the thermal reservoir which sets the 
base temperature T of the experiment (hereafter, the metallic electrodes), the heat has 
to travel across the sample of thickness e. Thus the dissipated power leads to a small 
temperature increase, containing a 2uj component 6T2 ~ p2- From the first order estimate 
SP{t) ~ {dPi{t) / dT)6T{t) which involves a product of two terms oscillating at u and 2u, 
one sees that a spurious contribution to the measured P3 comes from 6T2(t), i.e. from a part 
of the heating phenomenon. 

We shall see below that, close to Tg, this heating contribution to P3 vanishes in the limit 
of extremely thin samples (e — > 0). This is an illustration of the 'spurious' nature of this 
heating contribution. We emphasize that the heating we study here must not be confused 
with the 'heterogeneous heating' assumed in a model put forward recently by Richert et 



al. 
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17l |. This 'heterogeneous heating model' gives a phenomenological description of 



the intrinsic nonlinear effects in supercooled liquids, and was shown to account for nonlinear 
experimental data at uj in Ref. jl6| . In this model, each dynamical heterogeneity has its own 
fictitious temperature 'on top' of the temperature of the phonon bath. At the present time, 
it is the only model which allows to calculate the nonlinear susceptibility near Tg and we 
intend to give a thorough study of its predictions at So; in a future paper. To emphasize the 
difference between this heterogeneous heating model and the study of the present paper, we 
call the heating contribution to the nonlinear susceptibility investigated here 'homogeneous 
heating contribution'. 

The paper is organised as follows: in section [III we derive the relations which allow to 
calculate the heating contribution to P3. In section [Till we give the corresponding results, 
and clearly show the temperature and frequency ranges where the heating contribution is 
negligible in our experiments. We eventually show that in the high frequency range, i.e. 
when ojt{T) ^ 1, the heating contribution is damped by the finite relaxation time t(T) of 
the dipoles, and we estimate the magnitude of this damping. 
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II. HEATING EFFECTS CALCULATIONS 



In section III Al we establish the hnk between the temperature increase and the homoge- 
neous heating contribution to P3. The detailed calculation of the temperature increase is 
postponed to section III B[ 

A. From 6T to the homogeneous heating contribution to X3. 

This subsection is divided into three parts: we first establish an upper bound of the 
heating contribution to P3. We then move to an estimate of the damping of this contribution 
at high frequencies, due to the finite relaxation time of the dipoles. Last, we show what is 
the most natural quantity to plot to compare the heating contribution to P3 to the results 



of Ref. 



12|. 



1. The heating contribution without damping 

Let us consider a sample made of a supercooled liquid excited by an oscillating field 
Ecos{ujt). For small enough E values, the resulting linear polarisation reads : 

-^^^ = X'l cos{ujt) + x'l sm{ujt), (1) 

where eo is the dielectric permitivity of vacuum. In general, x'l and x'l strongly depend on 
frequency / = oj/27t. We define the frequency fa which characterizes the relaxation at a 
given temperature T, as the frequency where x'l is maximum. This relaxation frequency fa 
is strongly T dependent, and is of the order of 1/r where t(T) is the average relaxation time 
of the dipoles of the supercooled liquid. 

As evoked above, the volumic density of dissipated electrical power p{t) contains a d.c. 
term and a term oscillating at 2u (see below Eq. ( !T3|) ). The resulting heat propagates 
towards the 'thermostat'. The resulting average sample temperature increase ST[t) can be 
written: 

ST{t) = STo + ST2 cos(2wt - 02), (2) 
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where the mean dc temperature increase STq is larger than or equal to the ac one ST2, thus 
at any time ST{t) > 0. 02 is a phase shift related to heat transport that will be given in 
section IIIBI As our measurements [l2| give the nonlinear dielectric response averaged over 
the sample volume, ST{t) in Eq. is the temperature increase averaged over the same 
volume. Using Eqs. ([I]),©, we thus obtain for the nonlinear part of the polarisation due to 
heating effects: 

P{t)-Piin{t) fdx'i.^.A . , (d}(!{ \ . 

— ' -oT[t) cos{ut) + -——dT[t) sm(co't). (3j 



eoE \dT ' 'J V dT 

This is an upper limit of the heating contribution to the nonlinear response, since we have 
assumed that ST[t) induces instantaneously a modification of the susceptibility. As already 
advocated in Ref. [l^, this is questionable, specially in what concerns the contribution of 
ST2(t) which should be damped because of the finite relaxation time r of the dipoles. This 
point is adressed in the next section [Tl A 21 From Eqs. ([2]) and (jS]), one gets : 



Pjt) - PiUt) 



6T0 + -6T2COs{(f)2, 



+ -5T2sm(02^ 
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dT 
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cos(3u;)f:) 
sin(3ci;t). 



(4) 



The four terms in the right hand side of Eq. (jl]) give the nonlinear response of the system 
due to homogeneous heating: it contains two terms oscillating at u that we shall disregard 
since they contribute to the nonlinear part of Pi. We shall only keep the two terms of Eq. (j4]) 
oscillating at 3a; to obtain the heating contribution P^ h to the third harmonics P3. As we 
are interested in the heating contribution X3,h to the nonlinear susceptibility we define, as 
in Eqs. (4)- (5) of Ref. : 



P3 hit) E"^ . , , E"^ .. 

= ^X3,/» cos(3wt) + — X3,;,sin(3wt). 



(5) 



As a result, x'z^h ^^"^ x'i,h should not depend on E: their expression is given by identifi- 
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cation with the two last terms of Eq. (jlj) and defines, throughout this work, what we call 
the overestimated heating nonlinear susceptibility because it is obtained by neglecting the 
damping evoked above. 



2. Damping of heating contribution: an estimate 

We now move to the problem evoked above, namely the fact that the finite relaxation 
time r of the dipoles which contribute to the dielectric susceptibility should damp the mod- 
ification of this susceptibility due to the oscillating 6T2{t), specially in the case ujt > 1. As 
a consequence, the heating contribution to X3,h should be multiplied by a complex factor 
R{ujt), with a modulus \R{u)t)\ which is expected to be lower than 1. For a precise calcula- 
tion of R{ujt), one should replace Eq. ([3]) by an equation accounting for the dynamics of the 
dipoles in the case of a thermal bath where the temperature has an oscillating component, 
which is of great complexity. For an estimate, we make two very simplifying assumptions: 

(i) We assume that the dipoles have a Debye dynamics with a given characteristic time 
r(T). This is a simplifying assumption in the sense that when T is close to Tg, it is well known 
that Xi{^) is "stretched" with respect to a simple Debye law. In fact the Debye dynamics 
holds only at much higher temperatures, where the molecular motions are independent of 
each other, which allows to describe the non inertial rotational Brownian motion by the 
Smoluchowski equation for the probability distribution function of the orientations of the 



dipoles in configuration space IJ, ll8| . After an ensemble averaging of this equation,^ one 
gets the well known Debye equation for the dynamics of the average polarisation P [14,] : 

dP 

T— + P = eoAxiE cos{ut), (6) 

where Axi = = 0) — Xi{^ ^ oo) is the part of the static linear susceptibility 

corresponding to the slow relaxation process we consider. 

(ii) We assume that the main effect of the temperature variation 6T2{t) is to modulate in 
time the value of r while leaving unchanged the (Debye) dynamics. This can be justified by 
the fact that the temperature oscillation modulates the viscosity rj, thus also the relaxation 
time r which is proportional to r] [19\. Considering the temperature variations of Eq. 
T{t) is now given by 
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r{t) = run + (^) STo + 5T,(t), (7) 

where tu^ is the value of r at zero field. In the following, 6t2 will denote the amplitude 
of the 2u modulation of r due to 6X2(1) and corresponding to the last term of Eq. ([7j). 
Of course, using Eq. ([7]) for r(t) assumes that ST2{t) instantaneously fully affects r. A 
thorough modelization of this problem could lead to a more involved expression where 6t2 
would be weaker than in the above expression. As we shall find that the third harmonics 
is proportionnal to 5r2/r;j„, see below Eq. ([9]), we are led to the conclusion that our new 
estimate should, again, be slightly overstimated. 
We now insert r(t) in Eq. ([6]) and set : 

P{t) = Piin cos{ujt - + 6Pi cos(wt - *i) + P3 cos(3wt - ^3) + ••• (8) 

where Pun, "^Hn, ^Pi, ^i, P3, ^3 are to be determined. As we are only interested in the onset 
of nonlinear effects, Pun oc is much larger than 6 Pi oc and than P3 oc E^. This allows 
to neglect higher order harmonics (denoted by the dots in Eq. ([8])) and to resolve Eq. ([6]) by 
identification of the terms which have the same frequency and the same power of E. This 
yields: 



P3 
^3 



eoAxiE 6t2 



2 Tun ^l + ^UJTUnY^l + CiUJTUn? 
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^2 + arctan (urun) + arctan (3u;r/j„) + -, 



(9) 



where we remind that is the amplitude of the 2a; modulation of r due to 6T2{t) in the 
last term of Eq. ([7]): thus 6t2 oc E^, see Eq. ( IT^ below, which yields the expected P3 oc E^. 

We now have to compare with the result obtained if we start from Eq. ([3]) and use a 
Debye linear susceptibility. A straightforward calculation shows that in that case the third 
harmonics of the polarisation is given by the solution of Eqs. ([MH]) divided by the function 
R{ut) introduced above. The expression for R{ujt) is written using the complex notation : 



R[uT) = ^ — — X exp 
a/1 + (3wr;i„)2 



i X arctan 



-2uJTi 



lin 



1 + 'iiujTii. 



(10) 
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As expected R{ojt ^ 1) = 1 (no damping) and the damping arises at high frequencies 
since we draw from Eq. ffTOl) \R{ujt ^ 1)| < 1. More precisely \R{ut — )• oo)| = 1/3, which 
comes from the fact that r(t) enters in Eq. (jH]) as a factor of dP/dt: this gives, in Eq. 
(E]), a weight 3ut to P3, contrary to the case where one starts from Eq. where this 
weight is simply lot. Let us note that when the similar analysis is made for 6 Pi, by using 
the second term (oc 6Tq) of the right hand side of Eq. ([7]), no reduction is found at any 
frequency: The solution found for 6Pi{t) starting with Eq. ([3]), is exactly the same, in 
modulus and phase, as the one found by using Eqs. (jM?]). This shows that the reduction 
of the effect of ST2{t) on the polarisation comes from the fact that, in Eq. (E]), r and P 
oscillate together in time. Thus, Eq. f llOp can be seen as a first estimate of the fact that the 
dipoles damp the temperature oscillations, this damping being strong at high frequencies, as 
physically expected. Of course, one could build a much more thorough model of this effect, 
but the reduction given by Eq. (fTOj) will be shown to be quite realistic with respect to our 
experimental data (see below section UlI D[) . In practice, to compute the damped heating 



contribution to xs? "we first compute the overestimated contribution defined in section Hi A II 
and then multiply by the complex factor R{cot) defined above in Eq. ( ITOl) . 

3. How to single out the anomalous part of the nonlinear response 

Before moving to the calculation of the temperature increase, let us remind the re- 
lation that Bouchaud and Biroli predict [10], on quite general grounds, between xs 
NcorriT) -where Ncorr{T) denotes the T-dependent average number of djTiamically corre- 



lated molecules-. In Ref. 



10| . one finds the following scaling form for X3 '■ 



X3(C^,T) ^ '-^^^^^^N,orr{T)n{uJT) , (ll) 
Kb-L 

where the volume occupied by one molecule, and "H a certain complex scaling function 
that reaches its maximum at ut ~ 1 and goes to zero both for small and large arguments. 
This 'humped' shape of ["HI is due to the glassy correlations: In the 'no correlation case' 

HQ 

, NcorriT) Ti (ut) in Eq. ([11]) should be replaced by a function which reaches its 
maximum value at a; = 0. Thus Xsi^^ T) can always be considered as the product of a gen- 
eral prefactor e^i AxiY / ksT times a dimensionless term which summarizes the physics 
of the system. This is why a natural way to express the various contributions X3,i indi- 
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cates the kind of contribution) to xs is to divide them by this pref actor. We thus define 
a normahzed nonhnear susceptibihty Xs^j = / {eo{Axiyci^)x3,i- We shall consider the 
normalized heating contribution X^ h by dividing the heating contributions X3,h by the pref- 
actor. Clearly, when no heating contribution or any other spurious contribution is present, 
we expect X3(u;,T) = NcorriT)'H{ujT). 



B. Calculation of the temperature increase 



In this section, we now calculate the expression of 5T{t) that one has to introduce in 
Eq. (jl]) in order to obtain the nonlinear response in Eq. ([5]). 

The supercooled liquid is characterized by its thermal conductivity Kth and its specific 



heat c. As in Ref. 



20|, we consider that c is frequency dependent due to the fact that the 



20 



slow degrees of freedom cannot contribute to c for 
simplicity we neglect the small imaginary part of c 
depends neither on the frequency nor on the temperature T 
the plane of our copper electrodes [12 



requencies much larger than f^- For 



22|, and we consider also that Kth 



20 



22|. Let us define {x,y) as 
with 2; = for the lower electrode and z = e for 



the upper one. Due to their very high thermal conductivity and to their large thickness (6 
mm), the two electrodes can be considered, to a very good approximation, as a thermostat 



2l|. The temperature increase 69{x,y, z,t) of the supercooled liquid at point {x,y,z) and 
time t thus vanishes for z = and z = e. As the diameter D = 2 cm of the electrodes is 
typically one thousand times larger than e ~ 20 — 40/im (see we may consider that 66 
does not depend on {x,y). We obtain 69{z,t) by solving the heat propagation equation: 



dS9(z,t) 
' dt 



t^th 



d''se{z,t) 

dz"^ 



Pit), 



(12) 



where the dissipated power is given by 



p[t) = ^eoXiOjE"^ (1 + cos(2a;t - 0)) with (p 



— TT + 2 arctan 



X'l 



x'l - X'ii(^ 00) 



(13) 



This expression of p{t), where the fast, non relevant, degrees of freedom contributing 
to Xii^ ~^ 00) are separated from the slow degrees of freedom corresponding to glassy 
dynamics, deserves some comments. While the prefactor of the right hand side of Eq. ffT3|) 
is 'textbook' knowledge, the expression of (p is far less obvious (see the Appendix of |l5|). 
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This phase comes from the fact that the dissipation arises due to the friction of the dipoles 
with the surrounding molecules, and that this friction force T is proportional to dP/dt. 
The simplest example is the case of Debye dynamics where T is proportionnal to the first 



term of Eq. ([H]), see 



24| . Then, the power p corresponding to J-" must be given by J^v 



where v ~ dP/dt is the 'speed' of the dipoles. It follows that p ~ (dP/dt)"^ and by using 
Eq. ([1]), with the appropriate prefactor for p, the expression given in Eq. (fT3!) comes out. 
Note that if we consider the limit Xi ^ x'l ~ Xii^ ~^ oo); we get — )■ 0. This is the case 
of a metal where the response is not due to dipoles, but to electrons motion, and for which 
it is well known that there is no dephasing between p and E"^. On the other hand, in the 
limit Xi ^ Xi ^ X'li^ ^ oo), which happens in liquids when / ^ fa, one finds that the 2(jJ 
component of p and of E"^ are in phase opposition. In the case of a Debye dynamics, Eq. f lT^ 
can be rigourously derived |15)] because the expression of the friction force is explicit. When 
the dynamics does not follow the Debye 's law, which is the case of supercooled liquids close 
to Tg, the dynamical equation obeyed by each dynamical heterogeneity is not known jsj. 
However, Eq. ( IT3l) should still remain valid since the assumption J-" ~ dP/dt amounts to 
the lowest order development of the general idea that there is no friction if the dipoles do 
not move. 

Coming back to the heat propagation equation, due to the boundary conditions 69 {z = 
0,t) = = 69{z = e,t), Eq. f lT2|) is solved by decomposition in a series of spatial modes 
labelled by their wave vector K = mvr / e with m an odd integer. The mode m = 1 dominates 
the temperature increase, and we keep only the m = 1 and m = 3 modes since it is enough 
to get an accuracy of the order of 1%. By averaging spatially these two modes, we obtain 
the 6T{t) to be used in Eq. dH): 



Sn = 5T^{1 + ^] with 6T^ - 



34 y 24k 



th 



' Vl + (2a;r,,)2 3V1 + (2a;r,,/32)^ " 

where Tth = ce^/(Kt/i7r^) and 02,a = + arctan(2ci;rt/i) are involved in the dominant m = 1 
mode, while 02,6 = + arctan(2u;rt/i/3^) appears in the much less important m = 3 spatial 
mode evoked above 1231. 
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Following Eqs. ©-(E]), the two terms arising in ST2(t) in Eq. ( !T^ give a contribution to 
X3,h- These contributions are added, yielding the overestimated value of X3,h, as well as the 
damped value of X3,h (after multiplying by the function R{ujt) given in Eq. (ITU]) ). Then, as 
explained above, Eq. fITT]) is used to convert these values in terms of a contribution to X3. 



III. RESULTS 



A. Behavior at low and high temperature 



We shall first get some insight into the heating contribution X^^h to the total nonlinear 
normalized susceptibility X3 by extracting from the previous equations its frequency and 
thickness dependences. From Eqs. ((5]), (fT^ . we keep the leading term to obtain: 



IX. 



3,/i| 



dT 



(15) 



where we have defined Uth = ^nfth = l/(2Tth)- In Eq. f|T5|) two characteristic frequencies 
appear: Z^, which strongly depends on T, and the thermal frequency which is inversely 
proportional to the thickness of the sample. For a given experiment, and thus a given e, we 
have to distinguish two regimes: the low temperature regime where fa < fth and the high 
temperature regime where fa > fth- Besides, for glycerol, which is the liquid of interest 
here, the frequency dependence of Xi is Xi ~ f / fa below fa and x'i ~ if / fa)^^'^^ above /„. 
Last, \dxi/dT\ ~ {f/fX^ below /„ and \dxi/dT\ ~ {f/fay^'^ above /„. This allows to 
draw from Eq. f|T5|) the frequency and thickness dependencies of IX3 /1I. 

(i) In the low temperature regime, one gets : 





~ f-'e^ when < / < 






~ /-0.15g2 ^^^^ fa<f< fth 




^3,/i 


~ /-i-i5e° when < /. 


(16) 



In Eq. (fT6l) . the maximum over frequency of {X^^hl arises when the exponent of / changes 
its sign: this maximum value is thus proportionnal to e^, i.e. the heating contribution 
vanishes in the limit of very thin samples. This shows the non intrinsic (or spurious) nature 
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FIG. 1: (Color on line) Modulus of the overestimated heating contribution X3 /j for T = 204. 5K 
(triangles) and T = 217. 8K (circles), for glycerol and e = 19//m, as a function of the frequency 
normalized to the relaxation frequency /q(T). The open symbols are for a 'one sample experiment': 
the very flat maximum of |X3,h| at 204K is typical of the 'low temperature regime' (see text), while 
the 'sharp' maximum close to f = fa at 217. 8K characterizes the 'high temperature regime' (see 



text). The filled symbols correspond to the 'two samples bridge' described in Refs. ll|, [ij], where 
the heating contribution of the two samples, Cthin = 19/^m and Cthick = 41/im, cancel each other 
at low enough frequency, i.e. for / <C fthi^tMck) = 65Hz, see text. The net heating contribution 
is thus reduced in this setup: at 204. 5K this reduction is so strong that the maximum is shifted 
up slightly above 65Hz (the up arrow corresponds to 65Hz at 204. 5K) . At 217. 8K the weaker 
reduction is present only below 65Hz (see the down arrow). 

of the heating contribution which is studied in this paper. We note that for e ~ 19/im, 



which corresponds to the experiment of Ref. 12|, one typically gets fth = 300Hz. Since the 



standard definition of the glass transition temperature Tg corresponds to fa(Tg) = O.OlHz, 
the glass transition temperature is deeply in the low temperature regime. Therefore for 
e ~ 19/im, the heating contribution disappears close enough to Tg. This would not be true 
for samples where the thickness e is millimetric. 
(ii) In the high temperature regime, one gets : 



12 



[Deg: 


500 J 
400-' 


c 




TO 
Oi 


300-' 






CO 

v 


200- 






> 
o 

X"' 


100- 


o 


0- 


Phase 


-100 - 
10 



T=217,8K 
T=204,5K 




10-^ 



10" 



f/f 



10^ 



FIG. 2: (Color on line) Phase of the overestimated heating contribution X^ h for T = 204. 5K 
(triangles) and T = 217. 8K (circles), for glycerol and e = 19/xm, as a function of the frequency 
normalized to the relaxation frequency fa{T). The open symbols are for a 'one s amp le experiment'. 
The filled symbols correspond to the 'two samples bridge' technique of Refs. [ll|, L 



12l |. where the 



heating contribution of the two samples, Cthin = 19/im and ethick = 41/im, cancel each other at 
low enough frequency, i.e. for / <C fthi^thick) = 65Hz, see text. The net heating contribution in 
this setup thus depends on f / fth{ethick) ■ At 204.5K, / ~ fth{ethick) corresponds to ///„ ~ 13, 
see the up arrow, hence the separation between the filled and the open triangles ends only for the 
highest frequencies. Accordingly, at 217. 8K, the corresponding separation arises at f / fa — 0.3, 
not far from / ~ fthi^thick) which corresponds to f /fa — 0.13, see the down arrow. In the limit 
/ fthi^thick)-, the shift of the phase in the two samples setup with respect to the one sample 
setup is 7r/2. 



\X^A ~ /'-'e^ when Q < f < fth 
1^3,^1 ~ Z'-'e" when fth<f< fa 

\X^A ~ /"'-^'eO when < /. (17) 

Eq. ( IT71) shows that the maximum over frequency of \X^A arises for f = fa and is 
proportional to e°, i.e. independent of the thickness of the sample, just as the intrinsic non 
linear response. With fth ~ 300Hz, this high temperature regime onsets at T > 216K, nearly 
lOK below the maximum temperature studied in Ref. [l2|]. We thus expect that \X^A could 
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play a role for the highest temperatures reported in 12| . The open circles of Fig. [T] show 
the 'overestimated value' of {X^^hl for T = 217.8i^'. The exponents predicted in Eq. f|T7j) . 
for / < fth and f > fa, are well observed. Besides, by taking into account the prefactor 
not explicitly written in Eq. f[T7|) . one can check that the maximum over frequency of {X^^^l 
should arise for / ~ and give a value of order 0.15, as observed on Fig.[Il The comparison 
to the experimental data will be presented in section (jlll Cp . 



B. Heating contribution cancellation at low frequency with the bridge technique 

Before moving to the detailed study of X.sj, (a;, T), we investigate here the consequences 
of the fact that our experiment reported in Ref. 12| was performed with a bridge technique. 
Two samples of different thicknesses, ethm — 19/im and ethick — 41/im were used in a bridge 
to suppress the Su voltage due to the voltage source imperfection and to the (small) non 



linearity of the voltage detector 
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12|. Once the bridge is equilibrated, the voltages 
Vappiied,iuj applied to each of the two samples are stricly equal. The field applied onto the 
thin sample is thus larger than that applied onto the thick sample by a factor ethick/ ^thin- 
The subtraction operated by the bridge does not cancel the sought nonlinear response, since 
the latter goes as E^. We show now that this 'two samples technique' strongly reduces the 
values of {X^^hl as long as / ^ fth- Indeed, from Eq. f fT4l) . one finds, in the limit f <^ fth'- 

\5T2\ = ~ {VappHed,liuf , (18) 

i.e. the value of 6T2 is the same for the thin and the thick sample. By using Eqs. (I4]),([5]), 
this implies that the two heating contributions in the bridge setup perfectly cancel each 
other in the limit f <^ fth- 

This is illustrated in Fig. [1] where the filled symbols are the result corresponding to the 
bridge technique, while the open symbols are for a 'one sample' experiment. At T = 204K, 
one sees that \X^^h\ is strongly reduced by the bridge technique as long as / < 13/q, which 
amounts to / < fthi^thick) — 65Hz : the small remaining {X^^hl comes from corrections 
which depend on f/fth to Eq. flTB]) . and Fig. [2] reveals that the leading correction, in the 
limit / ^ fth, mainly produces a 7r/2 shift of the phase of X^^h- On the contrary, for 
/ ^ fth, one sees on Figs. [Tll2]that the heating contributions are similar, both in phase and 
magnitude, for a bridge setup and a 'one sample' experiment. 
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The same happens at T = 218K where / < fthi^thick) = 65Hz corresponds again to the 
frequency range where the bridge technique reduces (and phase shifts) the heating contribu- 
tion. Since the distinction between the 'high' and 'low' temperature regimes introduced in 
section fllll Ap involves a comparison between fth and /„, one concludes that the reduction 
of the heating contribution will be very important in the low temperature regime and much 
less important in the high temperature regime. More precisely, in a one sample experiment, 
for the low temperature range, the heating contribution has a very flat maximum around a 
few times fa- The example of the 204K curve of Fig. [1] shows that the two samples setup 
reduces the heating contribution so strongly that its maximum over frequency is shifted 
slightly above fthi^thick) = 65Hz, i.e. at / ^ fa- For the same reasons, the bridge technique 
extends slightly upwards the low temperature regime: the high temperature regime, charac- 
terised by a well defined maximum of X^^h located around / = fa, only arises at T > 218K 
instead of the 216K value evoked in the section (IIII Ap where the two samples technique was 
not taken into account. 

In the following, all the results shown in Figs. IMTOl correspond to a 'two samples' setup, 
in order to be directly comparable to our results of Ref. |l2[ |. 



3e compared to the experimental 
1^ . Two properties reported in 



C. Main features of the heating contribution to xs 

We shall now discuss the main features of the overestimated heating contribution 
(Figs. as well as those of the damped heating contribution (Figs. [SHS]). They will 

Xs (~ Ncorr'H) values at T = 210. 3K presented in Ref. 
are of interest here. First, when the temperature T is 
varied between 225. 3K and 194. OK, the maximum value maXfdX^l) over frequency of the 
measured IX3I increases by a factor ~ 1.5. Second, in this temperature interval, the fre- 
quency dependence of IX3I at each temperature, once rescaled vertically by its T-dependent 
max/ ( 1) value, fall onto a master curve depending only of f / fa- Besides, the phase of 
X3 also depends only of ///». This 'Time Temperature Superposition' (TTS) property of 
X3 is analogous to the similar properties of other well studied observables, e.g. Xi? for many 
supercooled liquids. Considering the fact that the above mentioned 1.5 factor is modest 
with respect to the nearly five orders of magnitudes displayed in Figs. we are led to the 
conclusion that the T variations of IX3I are hardly visible in these two plots. 
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FIG. 3: (Color on line) Modulus of the overestimated heating contribution l-'^s^/il for the two 
glycerol samples setup ethin = 19;um, ethick = 41/xm of Ref. [ij], as a function of the frequency 
normalized to the relaxation frequency fa{T). For comparison the measured values of IX3I [ij] are 
given for T = 210. 3K (filled red circles) . Several features of this overstimated heating contribution 
are at odds with those measured in [1^, e.g. the maximum over frequency increases with T, it 
occurs at a temperature dependent value oi f /fa and TTS is not obeyed (see text). 



From Figs. |3]|61 one gets the five following features, common both to the overstimated 
and damped heating contributions: 

(i) For a given fjfa, ^3,h mainly increases with T, at odds with the experimental behavior 
reported in Ref. [l2|. Besides the magnitude of the T-dependence of is much larger 
than that of X3 reported in Ref. 12|. 

(ii) The results for the modulus and the phase show that the heating contributions do 
not obey TTS. For the reasons explained in section IIIIBt the maximum over frequency of 
Xs^h arises for / ^ at low temperature, and for f c:^ fa in the 218K— 225K interval. We 
remind here that the maximum value over frequency of the experimental IX3I ^ \Ncorr'H\ 
arises at a very different frequency / ^ 0.21/^ -see Ref. 12[-. Summarising points (i)-(u), 
the frequency and T dependence of the heating contributions do not, at all, look like those 
of X3. 

(Hi) For / ^ 0.21/a, the heating contribution is always negligible with respect to the 



values of X3 reported in Ref. 12|]. More precisely, the heating correction to maxfdXsl) is 
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FIG. 4: (Color on line) Phase of the overestimated heating contribution X^^h for the two glycerol 



samples setup ethin = 19/^m, ethick = 41/im of reference [Ij], as a function of the frequency 
normalized to the relaxation frequency /q(T). For comparison the measured values of X3 at 
T = 210. 3K I12II are given (filled red circles): they do not significantly vary with T since the data 
reported in |12l | obey TTS (see text), contrarily to the heating contribution phases. 

smaller than 1.5% for T < 217.8/^, and reaches 5% for the highest temperature of 225. 3K. 
In this specific case of the highest temperature, subtracting the 5% heating contribution 
would decrease max/ (IX3I) by an amount equal to the error bar given in the Figure 3 of 
Ref.Q. 

(iv) For / ^ fa, the heating contributions decrease significantly faster with frequency 
(as f~^'^^) than IX3I. Thus the influence of the heating contributions disappears in the limit 

/> /a- 

(v) For T < 202 — 204K, the heating correction is always negligible for any frequency, i.e. 
a worst case analysis shows that it would amount to modify IX3I by less than the 5% error 
bar reported in Ref. This means that in the range [194K, 204K], the data of Ref. 
are absolutely free of any heating contribution, whatever the frequency. 
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FIG. 5: (Color on line) Modulus of the damped heating contribution X^^h for the two glycerol 
samples setup ethin = 19/im, ethick = 41/im of Ref. [12], as a function of the frequency normalized 
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to the relaxation frequency fa{T). For comparison the measured values of X3 at T = 210. 3K 
are given (filled red circles). The comments made in the caption of Fig. [3] still apply in this case, 
with the noticeable difference that the values of |^3,/i| are now smaller than those of Fig. [3] due to 
the fact that one takes into account the damping effect arising from the finite relaxation time of 
the dipoles (see text section Hi A 21 and Eq. (fTO]l ) . 



D. Evidence of a damping of the heating contribution to X'i 

In the range [194K,204K], the absence of heating contribution guarantees that Ncorr'H, 
obeying TTS, is well measured. In the [204K, 225K] range, the (Ncorr'H) giyc values of glyc- 
erol can be investigated by subtracting the heating contribution to the measured X^^fjT) 



reported at various temperatures in [12| . As we have seen, the heating contribution is weak 
or negligible in many cases, but it is interesting to study to what extent the expected TTS 
is better verified after such a subtraction. We define : 

(Ncorr'H) glyc = X3 — X^fi. (19) 

We performed this (complex) subtraction for the three temperatures above 204K in Ref. 



12| . i.e. 210. 3K, 217. 8K, 225. 3K. For each of these three temperatures, the frequency 



dependence is first fitted by a smooth function interpolating between the frequencies: these 
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FIG. 6: (Color on line) Phase of the overestimated heating contribution X^^h for the two glycerol 



samples setup ethin = IS^m, ethick = 41/im of reference [Ij], as a function of the frequency 
normalized to the relaxation frequency fa (T) . The comment of the caption of Fig. H] applies also 
here. 

three fits appear as continuous lines in Figs. [TllTOl Then, the subtraction of X3 /j to the 
measured X3 is performed, X3 being either the overestimated contribution(see Figs. [71 [8]) 
or the damped one (see Figs. [9l [10]) . The results of this subtraction appear as open symbols 
in Figs. [MOl 

When the overestimated value of X^^h is used, (Ncorr'H) giyc may differ strongly from the 
measured X3. Fig. [7] shows that the shape of the curve giving the modulus of (Ncorr'H) giyc 
vs frequency is modified at f /fa — 20 for T = 210. 3K. At 217.8K, a dip appears around 
f / fa — 2, and at 225. 3K the dip is present at f/fa — 0.8. These features come from the 
frequency range where the values of \X^^h\ are close to or larger than those of {X^l, as shown 
in Fig. [HI Depending on the difference between the phases of X3 and of X^^h, the effect on 
\Ncorr'H\giyc is more or less pronounced: at T = 217. 8i^ the dip at f / fa — 2 corresponds 
to a vanishing value of iNcorr'Hlgiyc, as shown by the vr jump on the phase (see Fig. [HI). 
At 210. 3K and 225. 3K, the non zero value of \Ncorr'H\giyc is related to a less pronounced 
effect on the phase (see Fig. [8]). To summarize, using the overstimated values of X^^h yields 
non TTS features in the resulting {Ncorr'H)giyc- We note that this is true already for the 
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FIG. 7: (Color on line) Solid lines: interpolated values of the experimental IX3I values reported in 



12l | for glycerol, obtained with the 'two samples bridge setup'. The intrinsic \NcorTT~i-\giyc of glycerol 
is plot as symbols and results from the complex subtraction — X^^h where the overestimated 
heating contribution of X3 (displayed in Figs. [3]- H]) is used. The resulting \NcorrT~l-\giyc curves, 
which display a dip or a strong variation of the / dependence for the three temperatures considered 
here, are non TTS, contrary to what is expected in the case of no heating effect. 

210. 3K curve corresponding to a temperature close to 204. 7K for which the curve of Ref. 
[1^ verifies TTS. The fact that our experimental X3(/, T) curves verify TTS is a strong 
indication of the overestimated character (already anticipated) of the heating contribution 
used to calculate the results presented in Figs. [T]- [HI We are thus led to the conclusion that 
the heating contribution has to be damped. 

We now perform the same analysis as above, using the damped X^ hif, T) values displayed 
in Figs. [5]- [6l At first glance. Figs. IMTO] reveal that for the damped heating estimate, the 
complex subtraction no longer produces the strong non TTS features depicted above. The 
resulting \Ncorr'H\giyc does not strongly differ from the original \X^\\ the value of the exponent 
giving the decay at ///a > 0.3 of \NcorrH\giyc is only 5% larger than the corresponding one 
for IX3I at 217. 87^ and 225. 3K. The above mentionned feature at 210. 3K and f / fa — 20 has 
almost disappeared. Last, the phase of (Ncorr'H) giyc is not significantly 'less TTS' than the 
original phase of X3. We thus conclude that the damped X3 values, introduced in section 
(III A2p and displayed in Figs [5]- [6l meet the requirement resulting from the last point of 
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FIG. 8: (Color on line) Solid lines: interpolated values of the phase of the experimental X3 reported 
in for glycerol, obtained with the 'two samples bridge setup'. Symbols: phases corresponding 
to the complex subtraction evoked in the caption of Fig. [7] (with the same symbols used). The vr 
jump occuring for 217.8K around ///a = 2 reveals that {Ncorry-lgiyc vanishes at this point. 



section ( jlll Cp . namely that (Ncorr'H) giyc should be TTS. This is an experimental indication 
that the damping factor introduced in Eq. f|TO|) is reasonnable, despite the two simplifying 
assumptions made along its derivation. 

We leave for future work the very difficult experimental task of isolating the heating 
contribution in itself. Two main ideas could be conceivable. First, one could try to push 
the 225K experiment in the rangef / fa ^ 1 where some non TTS feature might arise on 
the phase of (Ncorr'H) giyc if one extrapolates the calculations of Fig. [101 In practice it is 
extremely difficult to measure accurately X3 in this range of parameters, because of the d.c. 
heating contribution coming from the whole thermal circuit between each sample and the 
experimental cell 25|]. This d.c. heating contribution does not cancel in our two samples' 
bridge, and thus the balance condition changes at each value of the voltage source. The 
second idea would be to use the fact that the canceling effect of the heating contribution in 
the bridge arises only for / < fthi^thick)- As a result, the heating contribution to X3 should 
be different for two setups, say A and B, where ethick is different. In the frequency range 
fth{ethick,A) < f < fth{ethick,B), one may hope that the differences between the values of X^^h 
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FIG. 9: (Color on line) (Color on line) Solid lines: interpolated values of the experimental IX3I 



values reported in [1^ for glycerol, obtained with the 'two samples bridge setup'. The intrinsic 
\Ncorr'H\giyc of glycerol is plot as symbols and results from the complex subtraction X3 — X^^h 
where the damped heating contribution of X^ ^ (displayed in Figs. [5]- [6|) is used. The dips visible 
in Fig. [7]are no longer visible, and the resulting {Ncorr'Hlgiyc are approximately TTS. This suggests 
that the damped heating estimate is reasonable. 

could be detectable. Putting numbers, in the case ethick,A = 41/im and ethick,B = 25/im (see 
[isl). shows that this difference does not reach 10% of X3 (even for the optimal frequency), 
and is therefore very difficult to single out unambiguously. 



IV. CONCLUSION 



We have presented a thorough study of the 'homogeneous' heating contribution to the 
third harmonics experiments carried out on glycerol in Ref. [3| between Tg + 4K and 
Tg + 35K. We have emphasized the 'spurious' nature of this heating contribution by showing 
that it vanishes for thin enough samples and low enough temperatures. We have shown 
that the 'two samples' bridge technique, presented in Ref. [12|, widens the temperature 
interval over which the heating contribution is totally negligible with respect to the mea- 
sured Ncorr'H values reported in 121]. Globally, the heating contribution exhibits behaviors 
very different from those of the nonlinear normalized susceptibility X3 reported in \12\: It 
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FIG. 10: (Color on line) Solid lines: interpolated values of the phase of the experimental X3 



reported in [Ij] for glycerol, obtained with the 'two samples bridge setup'. Symbols: phases corre- 
sponding to the complex subtraction evoked in the caption of Fig. [9] (with the same symbols used) . 
The phase of the resulting {NcorrT~(-)giyc values are not significantly 'less TTS' than the original 
ones. This suggests that the damped heating estimate, displayed in Figs. [5]-[6l is reasonnable. 

(mainly) increases with T for a given ///a, its peak arises at much higher frequencies than 



m 



12| , its frequency dependence is faster, it does not obey TTS (except for the highest tem- 
peratures). At a quantitative level, we have shown that one can safely neglect the heating 
contribution to the maximum over frequency of IX3I for all the temperature interval studied 



m 



12l |. Thanks to the quite large temperature interval where the absence of any heating 
contribution is guaranteed, we have shown that the intrinsic Ncorr'H of glycerol obeys TTS 
in this temperature range. Extrapolating this TTS feature up to 22'bK allows to put an 
experimental constraint on the homogeneous heating contribution. We obtain the following 
important result: The fact that the homogeneous heating contribution must be damped 
because of the finite relaxation time of the dipoles is confirmed by an investigation of the 
shape and TTS property of the X^{f,T) curves. In addition, we may conclude that the 
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damping factor given by Eq. (fTOj) is reasonnable. 
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The difference lies in the fact that, because of the 'skin effect', the temperature oscillations 
are efficiently damped in the thick electrodes, while the d.c. component is not damped at all. 
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